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Abstract: We analyze multi-field inflationary systems which yield strongly scale 
dependent non-Gaussianity with a shape that is very close to the local shape. As 
in usual multi-fleld models, the non-Gaussianity arises from the non-linear transfer 
of scalar fleld fluctuations to curvature perturbations. Here we consider models 
in which higher order terms (loops) dominate over the lowest order source of non- 
linearity. The magnitude of non-Gaussianity depends on an infrared cutoff which is 
determined by our observational probes measuring non-Gaussianity. In our models, 
the running is positive and large {ni^c ~ 0.2) on CMB scales. The magnitude of the 
bispectrum is maximally of order 0(100), and grows on small scales. This can lead 
to interesting signals for large scale structure. 
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1. Introduction 

With the advent of precise cosmological data, it is now possible to constrain models of 
inflation by the measured magnitude and scale-dependence of correlated temperature 
perturbations in the cosmic microwave background (CMB) and from tracking density 
perturbations in dark matter from measuring the Large Scale Structure (LSS) of 
our universe. In these observations, it is found that the primordial perturbations 
coming from inflation are Gaussian to a remarkable accuracy, in agreement with the 
predictions of most single field models of inflation. 

Non-Gaussianity (NG) can be quantified by the magnitude of the bispectrum 
denoted /tvl (this is usually quoted at the equilateral point in momentum space 
where all three momenta are equal). For most slow-roll models, Jnl is smaller than 
1 [1, 2]. By comparison, the most recent constraints from WMAP5 [3] data are 
— 4 < fiyL < 80 for the local shape and —125 < /^^™ < 435 for the equilateral 
shape [4]. The Planck satellite is expected to improve the bounds to A/tvl < 7 [5]. 
There are also a large number of running and upcoming experiments probing LSS 
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scales (such as LSST, DES, SDSS, etc.) and they may allow us to eventually probe 
non-Gaussianity on smaller scales. 

In this note, we shall consider multi-field models with a large bispectrum (three- 
point correlation function) that is strongly scale dependent^. The running is positive 
(or blue which means that the NG grows as k increases) and can be achieved while 
keeping the power spectrum nearly scale invariant. It arises from loops (or higher 
order terms in the local ansatz) and the shape of the bispectrum is very well ap- 
proximated by the local shape multiplied by a logarithm. We provide a consistent 
setup where the 1-loop effect dominates the bispectrum while giving a sub dominant 
contribution to the power spectrum, and where higher loop contributions can be ne- 
glected. Since the running is positive, we can engineer a set-up where the curvature 
perturbation on CMB scales are extremely Gaussian while having a detectable NG 
on LSS scales. 

Running NG has already been considered in the context of DBI inflation [13, 14]. 
This model can have a strong NG signal due to a small and varying sound speed for 
the inflaton fluctuations [15]. The amplitude of the 3-pt can strongly run with scale 
if the sound speed varies but the running of the sound speed is exactly cancelled by 
the quickly varying Hubble constant along the trajectory. This is the key point of 
this type of model where the potential is steep but the inflaton moves slowly because 
of a speed limit. This causes the power spectrum to be scale invariant while the 
bispectrum can run wildly [16, 17]. 

The prospect of detecting large NG with large scale structure data has spurred 
much activity recently. LoVerde et al [18] have examined the possibility of using 
cluster counts and the galaxy bispectrum to constrain running /tvl. It was also 
realized in [19, 20], that NG of the local shape can induce a scale dependence of the 
galaxy/halo bias (see also [21, 22, 23, 24, 25, 26]). This effect can be easily found in 
the data and it results in a competitive bound on NG with local shape —29 < Jnl < 
70 [21]. At the time of this writing, there exists no signiflcant experimental bound on 
the running of NG with scale. Recently, Sefusatti et al [27] argued that Planck could 
bound riNG, the running of non-Gaussianity, with a precision Aung ~ 0.1 (0.3) for 
a local (equilateral) shape of non-Gaussianity. 

In our models, we flnd NG with a (nearly) local shape with a scale dependence 
such that the NG signal grows on small scales. The magnitude of the bispectrum 
grows with k with a model independent running of ung ~ 0.2 at CMB scale and 0.1 
on LSS scale. The strongest constraint on the magnitude of NG arises from n^. We 
flnd that Jnl ~ 100 can be achieved in principle. We also calculate the trispectrum 
tnl, which also runs. Before getting into the details, we summarize the basic idea 
and results. 

^ There has been much recent work in calculating the bispectrum and trispectrum in multi-field 
inflation, for some recent references see [6, 7, 8, 9, 10, 11, 12]. 
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2. Scale Dependence from Loops 

Local shape NG can be obtained in multi-field models of inflation, where each field is 
Gaussian but a non-linear relation between the inflaton perturbations and curvature 
perturbations induces NG. The original definition of the local ansatz for the curvature 
perturbation was done in real space [28] 

3 

C{x,t) = Ccauss + -^fNEiCcauss " {Ccauss)) > (2-1) 

where (causs is the Gaussian piece of the curvature perturbation, /tvl in this formula 
is by definition scale invariant. In momentum space, the above ansatz leads to the 
following bispectrum 

3 

= i2ny6%j:h)^^fMV^r^ , (2.2) 

where (C^O^a denotes a convolution, is the power spectrum (which is assumed to 

k^ 

be scale invariant, for simplicity) and defines the local shape. Many multi-field 

1 J- i 

models (such as curvatons [29, 30]) have local scale invariant NG of this type. The 
NG can also be scale dependent even if the shape is nearly local; for example, this 
is expected to happen when the NG is generated throughout the whole trajectory as 
opposed to simply at some fixed later time, such as in curvaton models. A particular 
model with this feature was considered by Byrnes et al [31, 32], where the scale- 
dependence arises from the dependence of /tvl on the (time-dependent) slow-roll 
and Hubble parameters. In their case, the NG decreases on small scales. 

We instead look for scale dependence coming from loops and higher order terms. 
Indeed, it was realized early on [33] that an additional contribution to the bispectrum 
in the ansatz Eq. (2.1) comes from 

(Ck^CkM,) = (^/A^^)'((CK)g,(CK)..,(CK).^3> • (2-3) 

This higher order contribution to the bispectrum has a structure similar from a loop 
contribution as it involves an integral over internal momenta. The integral converges 
in the UV but contains IR divergences if the power spectrum is nearly scale invariant. 
One can 'regulate' this divergence by introducing an IR cutoff in momenta 1/L ^. 

^These loops have been called c-loops [34]. They must not be confused with q-loops, or loops 
coming from the expansion of the quantum evolution operator prior to horizon crossing [35] . There 
has been much discussion recently on the physical significance of the IR divergences in loop calcula- 
tion in inflation. For c-loops, this IR cutoff is physical and depends on the observational probe and 
on how we measure the zero mode of curvature perturbations. We will justify this point of view in 
more detail in Sec. (3.1). 
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Doing so, the shape of this term is close to local up to a log [33, 36] 



(2.4) 
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If this term dominates the bispectrum, we will have a scale dependence with a running 
of order ting ~ As we will show later, the cutoff L is well approximated by the 
size of the universe today such that In/cL ~ 5 around CMB scale and ung ~ 0.2. 
The NG grows with scale becoming more important for smaller wavelength. Needless 
to say this is the interesting case as it gives rise to a stronger signal for LSS. 

Recently, Cogollo et al [37] and Rodriguez et al [38] have argued that loops can 
dominate in a particular 2-brid model. While their idea is very similar to what we 
propose, their particular model suffers from a problem pointed out in [31]. One of the 
fields that is assumed to follow a smooth classical trajectory is actually dominated 
by its quantum fluctuations, undermining part of their analysis. 

As we will show, the field that gives rise to NG in our model is also dominated by 
its quantum fluctuations. But this field plays no role in the inflationary trajectory and 
there is no inconsistency. We consider multi-field models of hybrid inflation where 
the inflationary trajectory is dictated by a single field but the surface of reheating 
(determined by when an extra waterfall/tachyon field starts condensing) fluctuates 
due to two fields [39, 40] (as originally envisioned by [41, 42] - see also [43, 44, 45] 
for similar models). 

In section 3, we describe the detailed set-up for the model, and describe the 
infra-red momentum cutoff. In section 4 we compute the power spectrum, and in 
section 5 we compute the bispectrum and trispectrum. We conclude in section 6 
with a discussion of these results. 



A simple way to move beyond single field slow-roll and generate NG is to have 
multiple fields. This type of model can quickly become very complicated and in 
order to simply illustrate the main physical effect of interest (namely large scale 
dependent NG from loops), we will consider a very simplified set-up. More general 
models and in-depth analysis of the model we present is left for future work. Consider 
a model of hybrid inflation with two real light scalar fields (0 and x) and a waterfall 
field T which ends inflation when it becomes tachyonic and condenses. In this paper, 
we will consider a rather general action, a more detailed and worked example is given 
in Appendix A. The action is (we follow the notation of [40]): 



3. Multi-Field Model 




V = V^nf(0) + Vudix) + Vr. 



mess 



(0,X,T)]. 



(3.1) 
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The only coupling between and x are through the tachyon which acts as a mediator 
or messenger. The form of Kness is taken to be 



K.ess(xTV(0,X) + C?(n 



n > 2 . 



(3.2) 



The function / interpolates from large and positive values (in Hubble units) during 
inflation to negative values after the system crosses a critical line in field space. 
Therefore during inflation, T has a large positive mass, its vev is driven to zero and 
its potential vanishes. Because of its large mass, this field will not fluctuate and it 
can be integrated out of the theory. In this model, inflation ends suddenly when the 
mass of the tachyon vanishes, which occurs on a line in fleld space parameterized by 



where the index "e" denotes the value of the flelds at the end of inflation. During the 
inflationary phase, and x have no direct coupling. To simplify further, we assume 
that Vhidix) ^ ^nf{(p) and we refer to as the inflaton from now on. The Hubble 
scale is then approximately given by 



and X is a "hidden" fleld during inflation which fluctuates but without much impact 
on the total energy density of the Universe. Nevertheless, its quantum fluctuations 
are still important as they will be felt as ripples on the surface of reheating. Indeed, 
at different point in space, the (slightly) different value of Xe will mean different 
critical value (pe for the inflaton resulting in more or less inflation in these different 
regions. This correlates directly in curvature perturbations (See Fig. (1)) Since the 
quantum perturbations of x mainly affect the surface of reheating, this system is 
well amenable to analysis through the 6N (or separable universe) formalism [46]. 
The idea is that the curvature perturbation on large scales is simply given by the 
perturbation in the number of efolds for each trajectories 



where the curvature perturbation ( is given by fluctuations of the scale factor a{x, t) = 
a{t)e'^^^'^^ and the difference in number of efolds is from a initial flat hypersurface 
to a uniform energy density flnal hypersurface. This formula does not take into ac- 
count possible interactions between the various flelds inside the horizon (on small 
scale) and it is only valid after horizon crossing where the evolution of the curvature 
perturbation is classical^. 

^The SN formalism will not account correctly for multi-field effects for modes inside the horizon. 
In our case, because the fields are uncoupled during infiation, we can solve for S(j) and Sx are horizon 
exit independently and follow the subsequent evolution of ( with the 5N formalism. 




(3.3) 



ax,t)=6N{x,t) , 



(3.5) 
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Figure 1: This figure depicts the trajectory in field space. The blue (dashed) line denote 
the surface of reheating defined by /(</)e) Xe) = and it is assumed to be thin. The classical 
trajectory is in the (j) direction (red/dotted line) but both 5(j) and 5x will induce curvature 
perturbations. 



The surface where inflation ends Eq. (3.3) is not a uniform energy density hy- 
persurface and a correction term must be included as discussed in [47, 44]. The 
correction term is very small in the hybrid scenario where the potential is very flat 
and it will be dropped in what follows. The number of efolds is given by dN = —Hdt. 
For the case where the classical trajectory is determined by a single field 0, one has 

N = - -jdcP' , (3.6) 

where the critical value of depends on the value of the field x the end of 
inflation (we dropped the subscript e and x = Xe unless otherwise specified^) and 
* refers to horizon crossing for a given mode. By varying 0* ^ 0=,, + 6(f) and then 
Mx + Sx) = <Pe + l6x + iJxV^ + ■■■ with 



7(x) 



dx 



(3.7) 



where we denote the zero mode of x by that is xi^^t) = x(^) + ^xi^^t) (for 
notational simplicity, the bar is omitted in any derivative subscript). We get at 
second order (using ^ = — A^') 



6N 



1 



N"^Hx' 



(31 



2 le 2 

where ' denotes derivatives with respect to 0. This can be reproduced using the 
formula of Vernizzi and Wands [47], for the case -C e*^ albeit they implicitly 

*The field x is evolving stochastically and the value of the field at the end of infiation is the sum 
of all fiuctuations created for each mode as they exit the horizon. 
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assume that all fields obey their equation of motion which is not true here for the 
field X- It is simple to show that A^' = dN/d(j) = 1/V2e^ Mp where the slow-roll 
parameters are 

= (^) \ e- = \mI ' . (3.9) 

The terms with N" involve derivatives of slow-roll parameters and will therefore be 
suppressed. To simplify the formula and the analysis we will consider the case where 
the slow-roll parameter at horizon crossing and at the end are equal, = . This 
is not true in many models and we will discuss at the end how that would affect our 
results. We thus drop all subscript referring to the time of evaluation. The mean of 
Eq. (3.8) is non-zero and as it is we will generate a one-pt function. To ensure that 
the mean is zero we can subtract a constant piece (keeping only the leading terms) 

C = - N'^bx - \n'iJx^ + i^x') , (3.10) 
which is of the form Eq. (2.1). 

C = Cl+C2-(C2) . (3.11) 
Note that this series terminates if 

1. the function 7 is such that '-f^y.^ and higher derivatives are small. 

2. A^" and higher derivative contributions are small. 

In this type of model, the function 7 could be anything and in the case where 
7,x^X > 7 the quadratic piece in 6x will dominate over the linear piece (in 6x) which 
ensures that the loop contribution to the bispectrum will dominate 

{C')^l'A(^xy) , (3.12) 

as we advocated earlier. In order for the power spectrum to be nearly scale invariant 
we will still need the 50 piece to be the dominant contribution to the power spectrum. 
There is no contradiction since the linear perturbation in does not contribute to the 
bispectrum (or gives a very small slow-roll suppressed contribution). Furthermore, in 
the case where the higher derivatives of 7 are suppressed, the higher loop contribution 
can be neglected, ensuring a consistent truncation. 

Another important point is that for the loop to dominate, the zero mode of x at 
the end of inflation [x^ which is the mean averaged over the size of the universe at 
the end of inflation) must be smaller then the l-cx deviation value of the perturbation 
around the mean. Taking the quantum perturbation to be of order 6x ~ H, we must 
have Xe < ^X- This is better seen in a specific model such as the one presented in 
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Appendix A. There, we use a model where 0e = fix^) such that 7 oc x and 7^^ ~ est 



and the series truncate. In such models it is clear that the quadratic term dominate 
over the linear piece when 

^-^^-.^^>1. (3.13) 



7 X 

It is then clear that the field x has to behave stochastically and is in no way following 
a classical equation of motion. The fact that x has essentially no effect on the 
inflationary dynamics prior to the reheating tells us that the stochastic behavior is 
unimportant during inflation. The value of Xey being stochastic, could have any value 
and it is therefore a free parameter. Before going into more details of the calculation, 
we need to discuss the choice of IR cutoff in the loop calculation. 

3.1 The IR cutoflF 

There has been much discussion in the literature about the choice of cutoff that 
should be used in loop calculations. For the calculation of quantum loops in the 
in-in formalism (prior to horizon exit), the correlations function of scalars appear 
to be sensitive to this choice of cutoff, and there is no clear understanding of how 
this cutoff should be set. But for the c-loops which we consider in this paper, the 
situation is considerably simpler and there is a natural choice of cutoff [48]. We will 
define the observed zero modes of the fields 0, x as 

I rL/2 ^ .L/2 

^o = j3 d^x(j) , xo = j^ d^xx, (3.14) 

where L is the largest scale over which we have measured the fields. The perturba- 
tions of the fields are then defined as 6(f) = (j) — (f>o, 6x = X ~ Xo- 

When computing correlators of 6N, we are actually interested in the correlations 
functions of the perturbations e.g. {6(f)^^5(f)j:^. From the definition of the perturba- 
tions, we see that the effect of subtracting the zero mode is to remove all Fourier 
modes with momentum k > L~^. Hence (50^ = 0^ for k > L~^, and zero otherwise. 
Similarly, we find 

The effect is to include a cutoff on any momentum integral. Due to the cutoff, 
the correlation functions will have an explicit dependence on L. This can be traced 
back directly to the fact that we are calculating correlation functions of perturbations 
like 6(j) = — 00, which have a direct dependence on L through (pQ. In this formalism, 
it is clear that all the dependence on L comes from the variation in the zero mode 
as a function of L as was discussed in more details in [48] (see also [49]). 
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To summarize, there is a natural cutoff L determined by the biggest scale on 
which we are able to measure the background zero mode of curvature. This is max- 
imally the size of the universe today L ~ 1/Hq. This coincides with the lowest k 
perturbations that are possible to observe now. Since there are about 5 efolds be- 
tween when the lowest observable wavenumber leaves the horizon and when CMB 
scales leave the horizon, we have kcMsL ~ e^. LSS are about two orders of magni- 
tude greater than CMB scales, giving k^ssL ~ e^°. 



4. The Power Spectrum 

We will first consider the two-point function (CfciCfc2)- Fo^^ the scalar fields, we have 

(5x|> = (50|> = (27rf 5=^(5^ ^O^fc), 

i 

m = %f . (4.1) 

and we consider a model where these expectation values are approximately constant 
and where V is scale invariant (independent of k). We will also assume that any 
intrinsic 3-pt functions are negligible, ((50^)") ~ and ((^Xfc)") ~ for n odd. In 
Fourier space the curvature perturbation is given by (from Eq. (3.10)) 

1 f d^P 1 

C.^ = N'5<j^^, - iV'75x.. - 2^'7,x J J^.^Xi.y^x^w + 2^'7,x (5x|> , (4.2) 

The "tree-level" contribution to the power spectrum arises from linear terms in the 
expansion of 5N , and it is easily seen to give 

= iV'2(l + i^)(2i) k,)P(k) . (4.3) 

i 

However, there is also a "one-loop" contribution which arises from the non-linear 
terms in the 5N expansion which leads to 

{'^ki^h) loop ^ ^ j (^^fci-fc'^^fc'^^fc2-fe"'^^fc") ' 



N-4(2nmy:k.) f p-m^ , (4.4) 



where the factor of (2) is from the combinatorics. For a scale invariant power spectra 
V, the integral is approximately 

^ d^k' 1 , , 

TTT^^ ^ > 4.5 

i/L \k - k'l^k'^ 
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where we use k as the upper hmit because for k' > k the denominator goes as fc'"^ 
with n > 3, and the integrand drops rapidly. The integrand has two simple poles 
which give logarithmic divergences. We regulate these by putting an IR cutoff on 
the integral. Hence for this example, we get 

^ d^k' 1 \n(kL) , , 

^ ^ 2—^ — - . (4.6) 

i/i (27r)3|A;_A;'|3fc'3 27r2 

This contribution will depend on the IR limit of the momentum integration. This 
limit is given by the size of the observable universe today, L ~ H^"^ as we discussed 
in Sec. (3.1). Modes of longer wavelength are already summed in the background 
value of the field. We thus find 

i 

Combining these terms yields 

(C,^> = (2-)^5^(E^^)^' (4-8) 

i 

= N'\2Trf5\Y, h)P [1 + 7' + l^V HkL)] . (4.9) 



We have defined the power spectrum for curvature with the superscript C,. The 

d\nk 



spectral index — 1 = ^7?^ is 



Note that the log contribution is positive (blue) and if this is the only contribution, 
we cannot match to the currently observed value of Ug ~ 0.96 [3]. For now, we simply 
impose that the log contribution contribute no more than a percent correction to Hg 

l^V < 10-2 , (4.11) 

which in turn implies that the non-linear contribution to the 2-point function must 
be subleading if log(/cL) ~ 1. 



5. Higher Point Functions 
5.1 Bispectrum 

We now compute the 3-point function {Ck^Ck^ChJ ■ Again, we find that this correlation 
function can easily be computed by expanding SN in terms of 6(j) and 6x- Since 6(f) 
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and 5x are Gaussian fields, the only non-trivial contributions will come from non- 
linearities in the 5N expansion. As in the case of the 2-point function, there is a 
natural separation into "tree-level" and "loop" contributions [50]. The contribution 
which is of lowest order in 7 is 

1 7 f d^k' 

= -^^N'-'{27Tf6^{yk,)^{27T^V)'^^ . (5.1) 
V 7 Ui kt 



The next term in the 7^;^ expansion is 

'73 j (27r) 



-^^N"i27rf6'iyh)-^ I 7^ ( ^ + 7perms 

^ ^ > '>S^^ J {2TiY\mk,+k'm,-k'\^ 



(27r)3 \kri\k^ + k'\^\k2-k'\ 



17.3 



= -^'N'\2Txf5\y k,) I^(2n'VfB{h, h) • (5. 
^—^ o 7'^ 

i 

Now the loop integral involves two different momenta 

Diagrammatically this is equivalent to a triangular loop of scalars (see Fig. (2)). We 
note that near the poles at k' = 0, k2, —ki, we get logarithmic divergences which are 
cut off by the IR scale L. This logarithmic dependence breaks scale invariance. So 
our shape S is a function of three variables which we choose to simply be the norm 
of all three vectors ^1,^2,^3. An estimate of the shape can be obtained by simply 
evaluating the integral around each poles, cutting off the momentum integration in 
the infrared at scale 1/L. So for example, the integrand 

d^P 1 
^ ^ (5.4) 



(27r)3^/3|^"; + ^/|3|^";_^/|3 



has a pole around k' = 0, and the integrand falls off rapidly when k' becomes of the 
same order as ki or k2. Hence we can approximate the integral around that pole as 

d^k' 1 ln(Min(ki,k2)L) 

■ + ■••. (5.5) 



(27r)3 + k'\3\g-k'\s 27i^kfki 



The same thing can be done for the other poles and for the various permutations. 
There are also points in parameter space where the integrand has a pole of order 




Figure 2: The 1-loop diagram. In our case, each vertex is accompanied by a factor of 
N'^^'^^ while each internal propagator is given by ^"^3^ . More detailed Feynman rules for 
use with the 6N expansion (which we are not carefully describing here) can be found in 
[51]. 

— * — # — * 

4. These poles occur in the squeezed limit where ki = — ^2 and hence 0. 
This shows that the bispectrum diverges in the squeezed limit, as is usual for the 
local shape. In principle, we can only measure k to a resolution ~ 1/L and the 
bispectrum, while large, is finite and of order L^/{3kf) in this limit. Hence the 
stronger poles that we have neglected are only important in the squeezed limit and 
they give contributions of the same order as the log terms in that limit. The full 
approximative shape is 



27r2 V k(k 



3j,3 
2 



The 1/3 term is only relevant for scales smaller than k ~ e ' j^- For larger k the 
shape is very well approximated by 

B{k,,k,,k,) ^ ^ln(Min(kOL)^ . (5.7) 

We show numerically in Appendix B that this is a good approximation. In Figure 
(3), we plotted the shape given by Eq. (5.6) in term of the usual variable X2 = k2/ki 
and X3 = k^/ki. When the bispectrum is scale invariant, ki is fixed to 1 (arbitrarily) 
but here we plotted the shape for different value of ki. As the figure clearly shows, 
the graph is very close to local and the magnitude grows as ki increases. At the 
equilateral point ki = k2 = k^ = k, the loop contribution to the bispectrum simplifies 
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Figure 3: Plot of the approximate shape B{ki, kiX2, kixs)x2x'^kf (with B{ki, k2, k^) given 
by Eq. (5.6)) in terms of X2 = ^ and X3 = ^ for ki = 0.5 (left) and ki = 1.5 (right). 
The shape was restricted to be in the quadrant defined by ki(l — X2) < kix^ < kiX2 due 
to momentum conservation and to avoid overcounting identical triangle configurations (see 
[52]). The shape is clearly very close to local with the strongest signal in the squeezed limit 
when /c3 = kix^ — > 0. The overall magnitude of NG increases with the wavenumber ki or 
as we consider smaller wavelengths. 



to 

(C^> = -f iV'3(27r)^5^(X]^,)^ln(fcL)(2vr^)^p3A . (5.8) 

If we compare the standard parameterization for local non-Gaussianities (Eqns. (2.1) 
and (2.2)) at the equilateral point to Eq. (5.1) and Eq. (5.8) and using the approxi- 
mation ^ A^'^P, we have 

/..--^^(l + ^ln(fcL)p) , (5.9) 

where the first term is the tree-level contribution, and the second term is the one-loop 
contribution. 

In the case of the two-point function, experimental bounds on the spectral index 
required the loop-contribution to be subleading. But there is no such requirement 
for the bispectrum. The loop contribution will dominate if 

^V\n{kL) > 1 . (5.10) 

In this limit we have 

5(7^ v)-2 

IfNil ~ - HkL) < 1001n(A;L) , (5.11) 

6 N'V^ 
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where we have utiUzed the bound T^'P < 10^^ and the normahzation V^^"^ ~ 
_/V'p^/2 10~^ from COBE data. We thus find, in this scenario, that one can 
easily generate local non-Gaussianity which is not ruled out by WMAP5 and can 
potentially be probed at Planck. Note that the magnitude of the non-Gaussianity 
increases logarithmically with momentum, suggesting that non-Gaussianity can have 
an important impact on the formation of structure at smaller scales. If we define the 
running of /jvl at the equilateral point 



riNG 



dink 

one gets in the loop dominated limit 



(5.12) 



" h^ ■ 

In the limit where non-linearities dominate, the running of /atl is thus independent 
of A^', 7 and 7,^. 

5.2 Trispectrum 

As in the case of the 3-point function, the only non-vanishing contributions will arise 
from the non-linear dependence of SN on 5x, so we can ignore 50 fiuctuations. To 
simplify notation, we define 

Ck = ^^xu + BJ ^s^xj;-n'Sxp - b (5x|> , (5.14) 

where A — —N'^ and B — —^N'^y^^. The last term ensures that we only keep the 
connected part of every diagrams. The tree level contribution (the term of lowest 
order in S) is 

(CkSk^ksCk,) = / (^^^ {^Xk,SXkJXk,-k'SXk>SXk^^k''^Xk'') + 5 perm , 

= {27r)HA''B^5^ (J2 ^) [^(^1 + h)P{ki)P{k2) + 11 perm] , 

-^A^B-{2n)H'{X.k)^. (5.15) 
where we have used that ~ N''^V and the shape is given by 

— * — * 

with the notation kij — \ki-\-kj\. The magnitude of the trispectrum is usually given 
by two numbers (tatl and Qni) corresponding to two distinct shapes: 



54 

TNLT{ki) + —gNL{P^{k2)P^{k^)P^{k^) + 3 perm) 
zo 



(5.17) 
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The lowest order contribution thus corresponds to g^qj^ — and Tatl — AA^B'^ /N'^ . 
The 1-loop contribution comes from the following term 

/d^p ■ ■ ■ (f' p^' 
^^^^12 ^ i^Xk.-k'SXk' ■ ■ ■ 5Xk,-kiv5Xkiv) , (5.18) 



(27r)3 k'3\k^ - k'\3\ki + k2- ^'|3|^3 + ^'|3 



h 5 perm 



The integral over momentum is difficult in general, so we will only estimate its value 
at the equilateral point \ki\ = k 

r>-,„P = 16(2.)»BV (J:*.) (2.W^^ (5.19) 



and thus 



4B^ 

^NL = j;^iA' + 4B'VHkL)) , 



QNL^O. (5.20) 

We see that the trispectrum is dominated by the non-linear contributions in largely 
the same regime as the bispectrum. Given the bound from — 1, the maximum 
value for tnl in this loop dominated regime is 

TNL ~ fj^ HkL) < 10^ HkL) . (5.21) 

Interestingly, the bound from WMAP5 on this parameter is \tnl\ < 10^ while Planck 
is expected to improve this bound up to \tnl\ < 560. 



6. Conclusions 

We have studied a simple class of models in which non-Gaussianity is dominantly 
produced by higher-order non-linearities in the transfer of fluctuations from the fun- 
damental scalars to the curvatiu^c. These higher-order non-linear order contributions 
are often referred to in the literature as "c-loops", and can dominate the lowest or- 
der "tree-level" contribution in the limit where ^Vln(kL) > 1, where 7 and 7^^ 
parameterize the non-linear transfer of fluctuations. In particular, J'nl ~ 100 can be 
achieved in these models. 

We have also found in these models that the magnitude of non-Gaussianity is 
scale dependent, with njvo ~ 0.2 at CMB scales and unq ~ 0.1 at LSS scale. 
Interestingly, the non-Gaussianity of the bispectrum is stronger at smaller scales. 
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where it can potentially be observed by large scale structure experiments. The shape 
of our NG signal is very nearly local. Moreover, this class of models yields a non- 
trivial trispectrum (parameterized by tml) that also runs. 

A number of open issues remain. In our model, we have assumed that the slow- 
roll parameter e is constant throughout inflation. This was necessary in order to 
have an observable effect from the end of inflation, but it requires tuning and it 
leads to a very flat power spectrum. It would be interesting to either relax this 
assumption in our scenarios or to look at a completely different set-up where the NG 
is not generated at the end of inflation. We expect that we can relax this assumption 
since we could have a case where f^L is very small on CMB scales but grows to be 
detectable on LSS scales. We note though that D-term inflation with a Coleman- 
Weinberg potential (as illustrated in Appendix A) has a natural regime with the 
required flat potential, eg ~ e/. From an effective fleld theory point of view (and 
from string theory models such as [40, 53]), the real tuning is in keeping all other 
allowed terms (such as a mass term for 0) subdominant to the Coleman- Weinberg 
potential. 

We have also assumed that the fundamental scalars (0 and x) ^^^6 Gaussian, and 
that all non-Gaussianity is induced by the non-linear transfer of 5x fluctuations to 
the curvature. Non-trivial NG can also arise from non standard kinetic terms, or 
a steep potential for x (which unlike the inflaton does not have to satisfy slow-roll 
conditions). Loop corrections then have a richer structure although the basic idea 
remains the same. Of particular interest are models like DBI inflation where the 
spectral index is nearly one and entropy modes being converted to curvature at the 
end of inflation can also be observable [54]. This scenario has been analyzed recently 
in [55] based on methods developed in [56, 57] (see also [58] [59]) and a mixture of 
equilateral and local NG has been found. It would be interesting to consider the 
regime where the loop dominate in this kind of models. 
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A. A Specific Model 



The discussion in the text is very general and the ultimate goal of having dominant 
loop contribution in the bispectrum inducing a large running may be achievable in a 
variety of ways. Here we we illustrate the necessary ingredients with a specific model 
(based on [39]). Take an inflationary potential 



^2 



2 

Vcw{x) = {x^ + if \n{x^ + 1) - Inx^ + (x^ - 1)^ ln{x^ - 1) - 4 In 2 



(A.l) 



where = ^ has mass dimension 2 and A and g are dimensionless couplings. 
The reader will recognize this as the Coleman- Weinberg potential. There is a regime 
in parameter space where the inflaton does not move very much with 

0^ - 0^ = § (A.2) 

and the slow-roll parameter is also nearly constant^ 

'^Mn2VM^2 AV(ln2)^M^^ 
7^ J -2 2^^ • (^-^^ 

The and x power spectrum are simply given by V = and they will remain 
approximately constant until the end of inflation if r/^ and rj^. are much smaller than 
1. Now consider a simple potential for x 

V1.id = i^W^ . (A.4) 

This potential drives x to but the field will fluctuate and acquire some stochastic 
value X which in general will be non-zero (although small). The tachyon potential is 
of the form Eq. (3.2) with the surface of reheating defined by 

= /(0e, Xe) = -g^i + \m + PXI ■ (A.5) 

We choose a model such that Xe = Xe + 0e- Since the function / is quadratic 

in both fields, the transfer function is simply 

/5Xe 



7 



and 7 V ~ - while 7 ~ 0. Note that 7 and 7.y can both be either sign depending 
on (3. The curvature power spectrum is (from Eq. (4.8)) 



i/2 / , 72 H-2 kL . 



5 



By integrating the EoM of motion of 0, in the hmit a; 1, one can check that 0* ^ 0e is a 
good approximation as long as jyl^ ^ 2\/2 hi 2Afe ^i^gj-g jy^ jg number of efolds between horizon 
crossing and the surface of reheating. 
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We are interested in the regime where 



1 > — 

(27r)2 



(Ai 



and where the power spectrum P^, f^i and t^l are well approximated by 

e 



Vc 



f. 



NL 



TNL 



IT 

'6"(21n2)2A2M6 ' 

5 In 2(3^H'^M^\ii kL 
67r2A2c/2(27r)2^2 

p'^M^H^ In kL 
87r4(27r)2A2^2^3 • 



(A.9) 



Note that none of these observables depend on the precise value of x at the end of 
inflation although in order for the loop contribution it must be that the average value 
of the zero mode of x at the end of inflation is smaller than H. We show a point in 
parameter space (see Table (1)) where all the conditions mentioned in this section 
are respected. 



9 


A 






10-3 


2 X 10"^ 


0.1 


2.98 X 10-6 



Table 1: A point in parameter space. A was first chosen and was solved for by matching 
to COBE data. The parameter g is then constrained such that (p < Mp by at least two order 
of magnitude. /? is a free parameters that determine the magnitude of NG. The stochastic 
value of X at the end of inflation in this model is less than H although none of the observables 
depend on its precise value. One can check that for this choice of parameters: X'^ipe ^ ^'x^^ 
0.03 and 72 ~ 0.001. Also for this choice of parameter the potential is very 



-y^H'^ In kL 
' (27r)^ 

flat with e 



~ ~ 



10-" while e2 



10 



-29 



(for z^2 10 2 giving r/^ ~ 10 





In kL 




Inl 


tnl 


Us - I 


nNG 


CMB scales 


~ 5 


2 X 10-9 


1.6 


1152 


0.004 


0.2 


LSS scales 


~ 10 


2 X 10-9 


62 


4.5 X 10^ 


0.004 


0.1 



Table 2: Predicted value for various parameters. This point was deliberately chosen to 
illustrate the possibility of having a non-observable level of NG at CMB scale but with a 
very detectable signal for LSS. 

As can be seen from Table (2), this simple model can lead to interesting observa- 
tional signatures whereas the CMB is very Gaussian, but significant NG appears for 
large scale structure. On the other hand, not everything is perfect since the spectral 
index is nearly one in tension with the most current WMAP5 data. This is direct 
consequence of working with a model where ~ . If this assumption is relaxed. 
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a running will be induced but the non-Gaussian signal coming from the end of in- 
flation will also be reduced. This can potentially be compensated by varying other 
parameters but this required more detailed analysis keeping track of the time of eval- 
uation for each quantity. We leave this for further work. Note also that the tension 
between ~ 1 and WMAP5 data can also be reduced if cosmic strings (which are 
generically produce in these hybrid models) contribute to the density perturbation 
spectrum [60] . The cosmic strings will add their own source of NG which will further 
constrain the model [61]. 



B. Numerical Evaluation of the Integral 

The integral 

B(h, h, h) = I 7^ I ^ ^^-^ ^ + 7 perms 1 (B.r 

can be evaluated numerically in Mathematica. There are three poles in the integrand, 
and the IR cutoff discussed in Sec. (3.1) is most easily implemented by setting the 
integrand to zero whenever k' is within 1/L of a pole. 

On general grounds, one expects that the integral is can be written as 

B{h. h) = HhL)F^{\h\, Ihl) + F,{\hl Ihl) + . . . , (B.2) 

i 

where the neglected terms contain powers of 1/kiL. For momenta relevant to CMB, 
these terms are very small. Note that both Fq and Fi are homogeneous Lorentz 
invariant functions of the external momenta with degree -6 (i.e., are scale- invariant), 
and thus are completely determined by the norms of any two of the momenta. 

As argued in Sec. (5.1), the leading term is dominated by the poles of the 
integrand, and we expect it to be of the form 

< (n,*i) 

Because Fq i are scale-invariant, we can numerically integrate the shape at various 
scales and fit the result to our ansatz in order to determine the magnitude of Fo,i. 

For example, we numerically integrated B at the equilateral limit ki = k for 
various values of k near the range kL ^ 150, and we fitted the results to the ansatz 
B{k,k,k) X {kLY = co\n{kL) + Ci. This fit yielded the expected Cq ~ where 
the second term ci was ~ 5% of the first term. Similarly, one can integrate B for 
several external momenta in the squeezed limit where kiL = ksmaiiL ~ 150 and 
k2L = k^L = khigL ~ 10000. Since the local shape should dominate in this limit, one 
would fit this to B{ksmaiu hig, hig) X {ksmaiiLf{kugLf ~ do \xi{ksmaiiL) + di- Again 
as expected, one finds do ^ where the non-logarithmic term di was ~ 8% of the 
logarithmic term. We conclude that the logarithm captures the leading behavior. 
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